SCIENTIFIC 

REPORTS 




OPEN 



SUBJECT AREAS: 

SURFACES, INTERFACES 
AND THIN FILMS 

TWO-DIMENSIONAL MATERIALS 



Received 
2 August 2013 

Accepted 
1 November 201 3 

Published 
18 November 2013 



Correspondence and 
requests for materials 
should be addressed to 
J.J.Z. (zhaojj@dlut. 

edu.cn] 



From Boron Cluster to Two-Dimensiona 
Boron Sheet on Cu(l 1 1 ) Surface: Growth 
Mechanism and Hole Formation 



Hongsheng Liu 1 , Junfeng Gao 12 &Jijun Zhao 1 



1 Key Laboratory of Materials Modification by Laser, Ion and Electron Beams (Dalian University of Technology], Ministry of Education, 
Dalian 1 16024, China, 2 Beijing Computational Science Research Center, Beijing 100084, China. 

As attractive analogue of graphene, boron monolayers have been theoretically predicted. However, due to 
electron deficiency of boron atom, synthesizing boron monolayer is very challenging in experiments. Using 
first-principles calculations, we explore stability and growth mechanism of various boron sheets on Cu( 111) 
substrate. The monotonic decrease of formation energy of boron cluster B N with increasing cluster size and 
low diffusion barrier for a single B atom on Cu(l 1 1) surface ensure continuous growth of two-dimensional 
(2D) boron cluster. During growth process, hexagonal holes can easily arise at the edge of a 2D triangular 
boron cluster and then diffuse entad. Hence, large-scale boron monolayer with mixed hexagonal-triangular 
geometry can be obtained via either depositing boron atoms directly on Cu(l 1 1) surface or soft landing of 
small planar BN clusters. Our theoretical predictions would stimulate further experiments of synthesizing 
boron sheets on metal substrates and thus enrich the variety of 2D monolayer materials. 

Since the first experimental isolation of graphene by Novoselov and Geim in 2004 1 , there have been 
tremendous interests in the two-dimensional (2D) monolayer materials of only one-atom thickness 2,3 , such 
as h-BN sheet and silicene. Among them, graphene and silicene constitute a special kind of semimetals with 
Dirac cones at the Fermi level 4 6 , whereas infinite fr-BN sheet is a wide-band-gap semiconductor 7 . Certainly, the 
diversity of physical properties of these 2D monolayers is crucial for their wide applications in future nanoscale 
materials and devices. 

On the periodic table, boron is the neighbor of carbon with three valence electrons in the 2 s and 2 p orbitals; 
thus it is able to form sp 2 hybridization like carbon. Intuitively, boron atoms are expected to form fullerene-like 
cage configurations and graphene-like monolayer sheet. In 2007, Yakobson's group 8 predicted a very stable B 80 
fullerene cage analogous to C 6 o buckyball. However, hollow cages are indeed not the ground state structures for 
the medium- sized B N clusters from N = 68, whereas core- shell configurations are more thermodynamically 
preferred based on ab initio global search by our group 5,10 and others 11,12 . The energetic unfavorability of these 
empty boron cages can be attributed to the electron deficiency of boron, which tends to adopt more compact 
structures. 

Stimulated by the B 80 fullerene cage, a 2D boron sheet (BS), namely "a-sheet", with appreciable stability was 
constructed 1315 . By incorporating periodic hexagonal holes in the triangular lattice to reach the balance between 
three-center (3c) and two-center (2c) bonds, this mixed hexagonal-triangular boron sheet (ht-BS) is more stable 
than previously proposed buckled triangular boron sheets (f-BS) 16 ~ ls , which can be explained by a chemical 
bonding picture that the hexagon holes serve as scavengers of extra electrons from the filled hexagons 20 . There 
have been some recent activities of searching more stable structure for 2D allotropes of boron, including the snub- 
sheet, g 1/8 -sheet, g 2 /i5-sheet, o^-sheet, pVsheet, struc- 1/8-sheet 21 " 24 . All these ht-BSs can be constructed by 
carving different patterns of hexagonal holes within the triangular sheet and described by a hexagon hole density 
T|, defined as the ratio of number of hexagon holes to the number of atoms in the original t-BS. 

These boron monolayers of abundant structures certainly bring new members to wonderful 2D atomic films 
and provide novel building blocks for nanoscale materials and devices. But monolayer BS has not been experi- 
mentally synthesized yet due to the electron deficiency of B atom. In previous experiments 25,26 , single- walled and 
multi-walled boron nanotubes have been observed, which hints the possible existence of boron monolayer sheet 
since a boron sheet can be viewed as unrolling a boron nanotube. With metal passivation to stabilize the sp 2 
hybridization, similar silicene monolayer has been successfully fabricated on Ag(lll) and Ir(lll) surfaces in 
recent experiments 6,27 30 , providing important clues for BS synthesis. On the other hand, small boron clusters in 
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vacuum were proved to adopt quasi-planar configurations 3137 , 
which may act as precursors for experimental synthesis of BS on 
metal surface via soft-landing of cluster beams 38 . 

As experimentalists have to spend lots of money and time by 
testing numerous metal substrates and various growth conditions, 
it would be instructive to investigate the synthesis of BSs on metal 
surfaces from the theoretical aspect. Recently, Zhang and co-work- 
ers 39 proposed boron monolayers can keep their structures with 
slight buckling on several metal surfaces based on theoretical 
exploration. Using first- principles simulations, Yakobson's group 
suggested that stable BS with a certain density of hole can be derived 
by the deposition of B atoms on Au and Ag surface or saturation of B- 
terminated MgB 2 surface in a B-rich environment 40 . 

Typical epitaxial growth of a boron sheet on metal surface can be 
divided into three stages: (i) a certain amount of B N clusters form on 
metal surface by aggregation of B atoms or soft landing of these B N 
clusters; (ii) B N clusters keep growing into pieces of BS sheet by 
attaching more and more isolated B atoms; (iii) further annealing 
to achieve stable monolayer structure and to remove defects. As ht- 
BS is more stable than f-BS, the final stage would be important for 
reaching the optimal density and pattern of holes in ht-BS. Although 
previous theoretical simulations 39,40 have made certain progresses, 
some fundamental mechanisms of BS growth still have to be clarified. 
For example, what is the most stable configurations for the B N clus- 
ters on a specific metal surface at the initial stage of BS growth? How 
do B atoms attach to the edges of a seed B N cluster? How do the holes 
form and assemble into patterned ?if-BS? If various BSs have grown 
on metal surface, how to identify and distinguish them using experi- 
mental means? 

To address these critical issues, here we investigate the equilibrium 
geometries and interaction mechanism of boron clusters and sheets 
on Cu(lll) surface using comprehensive first-principles calcula- 
tions. Similar to the free-standing sheets, /if-BSs on Cu(lll) surface 
energetically prevail against both f-BS and h-BS. Interestingly, 
the embryo of ht-BS structural motif emerges as early as B n cluster 
on Cu(lll) surface and competes with those from f-BS motif. 
Moreover, boron atoms tend to congregate on Cu(lll) surface with 
small diffusion barrier to form 2D clusters. During growth of these 
2D clusters with triangular networks, hexagonal hole can easily occur 
and migrate, which promotes the formation of ht-BS on Cu surface. 
Our theoretical results provide valuable guidance of how to achieve 
2D boron sheets in future experiments. 

Results 

Boron monolayer on Cu(lll) surface. The atomic structures for 
various boron monolayer sheets on Cu(lll) systems are illustrated 
in Figure 1 and Figure SI in the Supporting Information. After 
relaxation, the BSs with slight buckling retain their elementary 
geometries, implying possible existence of planar boron structures 
on Cu substrate. The amplitudes of the buckled height h, defined as 
the vertical distance between the highest and lowest boron atoms 
with regard to the Cu surface, were calculated for eight types of BS 




Figure 1 | Atomic structures of two representative boron monolayers on 
Cu(lll) surface, (a), (b): g 2 /i5-sheet; (c), (d): triangular sheet, (a), (c): top 
views; (b), (d): side views. 



considered and summarized in Table 1. For all /jf-BSs, the buckled 
heights are around 1 A, comparable to that for silicene on Ag(l 1 1) 41 . 
Meanwhile, the buckled height is 1.2 A for the triangular sheet, 
which is slightly larger than those for ht-BSs. In addition, we 
measured the average distance d between boron atoms and Cu 
substrates. As shown in Table 1, the average B-Cu distance ranges 
between 2.1 A and 2.2 A, which is comparable to that in the known 
B-Cu alloys (2.1 ~ 2.2 A) 42 and suggests moderate covalent bonding 
interaction between boron and copper atoms. The buckled 
geometries and various hole distributions of BSs on Cu(lll) 
surfaces may increase the difficulty to identify their structures. To 
help experimentalists distinguish different BS monolayers in the 
future, here we provided simulated scanning tunneling microscope 
(STM) images for typical BSs on Cu(lll) surface (see Figure S2 of 
Supporting Information). 

The formation energies for various BSs in vacuum and on Cu( 111) 
surface are summarized in Table 1. Compared to the standalone 
sheets, the formation energies of boron monolayers on Cu(l 1 1) sur- 
face are reduced by about 0.2 eV/atom, indicating that Cu substrate 
can stabilize BSs, similar to the Ag(l Unsupported silicene sheets 41 . 
On Cu(lll) surface, the most stable BS is the g 2 /i5-sheet with Ef orm = 
0.207 eV/atom, while the g 1/8 -sheet and oc-sheet are also rather stable 
with only 0.013 eV/atom and 0.019 eV/atom higher in energy. 
Indeed, these three isoenergetic ?zf-BSs, obeying the rule of 1/9 £ 
T| £ 2/1 5 22 , are the most favorable ones in vacuum 13,22 . Also note that 
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g, / 8 -sheet 
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oc-sheet 
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0.99 
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triangular sheet 
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0.338 
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1.21 
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the formation energy of the fully triangular sheet on Cu surface is 
0.338 eV/atom, higher than the most stable g 2 /i 5 -sheet by 0.131 eV/ 
atom. Such noticeable energy difference between /if-BS and f-BS 
would provide the driving force for structural transformation from 
fully triangular to mixed hexagonal-triangular sheets, as we will dis- 
cuss in the Section 3.3 below. Test calculations with DFT-D2 
method 43 for three typical boron sheets on Cu surface have also been 
performed and showed that inclusion of dispersion correction only 
lead to a systematical shift (about 0.023 eV/atom) of formation ener- 
gies without changing the energetic sequence of different sheets. 
Therefore, standard GGA-PBE method is sufficient for describing 
the relative stability of the present systems. 

Using the structures reported in our previous study 41 and the 
present theoretical scheme, the calculated formation energy of sili- 
cene on Ag(lll) surface is about 0.2 eV per silicon atom, which is 
comparable to that for various ht-BS on Cu(lll) surface. In recent 
years, silicene has been successfully synthesised on various metal 
substrates 6,27-30 . Thus, we anticipate that boron sheet may also be 
synthesized on metal surface with the method similar to that for 
silicene synthesis. In nature, however, there is no boron allotrope 
with layered structure analogy to the graphite. The bulk phases of 
boron solids (such as cx-boron and P-boron) are mostly based on 
basic unit of B 12 icosahedron 44 . To compare the mixed hexagonal- 
triangular boron sheets (?if-BSs) with a dense 2D arrangement of 
icosahedra, a monolayer composed of B 12 icosahedra on Cu(lll) 
surface was constructured (Figure S2). Using the same theoretical 
scheme, the formation energy of Cu-supported B 12 -based sheet is 
0.214 eV per boron atom, which is almost energy degeneracy with 
the most stable ht-BSs (see Table 1). This result clearly indicate that 
the boron thin film would eventually become bulk-like composed of 
B 12 icosahedra if enough boron atoms are deposited on Cu surface. 
Note that their atomic densities of boron on Cu(lll) surface are 
different, i.e., 0.56 A' 2 for B 12 icosahedron sheet and 0.35 A" 2 for 
fef-BSs respectively. Therefore, in order to achieve boron monolayer, 
the amount of boron atoms deposited on Cu substrate must be 
strictly controlled within a small value (about 0.35 atom per A 2 ). 

Moreover, using the structures reported in our previous study 41 
and the present theoretical scheme, the differences in formation 
energies of different silicene@Ag superstructures are about 0.03 eV/ 
atom, which are comparable to that for various mixed hexagonal- 
triangular boron sheets on Cu surface (about 0.02 ~ 0.03 eV/atom, 
see Table 1). In experiments, coexistence of various silicene@Ag 
superstructures has been observed 28,45 . Consequently, the final boron 
sheet synthesized on Cu surface might be a 2D polycrystalline as a 
mixture of these ht-BS domains, i.e., different patterns of hexagonal 
holes on the filled triangular lattice. 

The adhesion strength between ht-BSs and Cu( 111) surface can be 
directly measured by the E c defined by Eq.(6), which ranges between 
0.053 and 0.066 eV/A 2 (see Table 1). According to Bader charge 
analysis, each boron atom gains only about 0.04 electrons from Cu 
substrate. This can be attributed to the small electronegativity differ- 
ence of these two elements (2.04 for B and 1.90 for Cu according to 
Pauling's definition). Such moderate B-Cu interaction would be 
beneficial for stripping boron monolayer films from Cu substrate, 
which is a necessary step for future utilizations of the 2D boron 
materials. In addition, the lowest-energy configurations of ht-BSs 
on Cu(lll) surface, i.e., g 1/8 -sheet, g 2/15 -sheet, a-sheet, are also the 
most energetically favorable ones in vacuum. Therefore, once an as- 
prepared boron sheet is stripped from Cu substrate, it can stand alone 
and remain stable during the transfer process between different 
substrates. 

It is known that the graphene-like honeycomb lattice of boron is 
unstable in vacuum due to electron deficiency 46 . Since the Cu-to-B 
charge transfer is too small to remedy the electron deficiency of 
boron, the h-BS cannot survive even on Cu(lll) surface and will 
collapse into a disordered sheet with triangular domains and large 



voids (see Figure Sl(m, n) of Supporting Information). This result is 
in accordance with previous finding by Yakobson's group 40 but in 
contrary to the prediction by Zhang et al. 39 , who suggested that h-BS 
on some metal surfaces (Mg, Al, Ti, Au, Ag) is more favorable than 
both f-BS and ht-BS. This difference might be attributed to stronger 
BS-substrate interactions (about 0.6 —1.2 eV/atom) from Zhang's 
calculations. 

In real experiments, boron monolayer must be synthesized at 
finite temperatures, where thermal perturbation may introduce some 
structural defects and thus disturb the geometry integrity of 2D 
atomic sheet 41 . To examine the thermal stability of boron sheets on 
Cu surface, we preformed ab initio molecular dynamics (AIMD) 
simulations for a-sheet on Cu(lll) surface as a representative. 
Using VASP code, AIMD simulations were carried out at 500 K 
and 800 K with time step of 1 fs. At both temperatures, no topo- 
logical defect was ever generated during the entire simulation time of 
6.5 ps (Figure S3), confirming the high thermal stability of boron 
sheets on Cu(lll) surface. 

Structures of boron clusters on Cu(lll) surface. From the above 
discussions, one can see that various ht-BSs are all stable on Cu(lll) 
surface. To further elucidate the growth behavior and nucleation 
mechanism of boron sheets, it is thus necessary to explore the 
geometries and stabilities of small boron cluster on metal substrate, 
as we did before for graphene and silicene clusters on metal 
surfaces 41,47 49 . Several critical issues have to be addressed: (1) do 
small boron clusters prefer two-dimensional or three-dimensional 
structure on Cu surface? (2) what is the structural difference for 
boron clusters in vacuum and on metal substrate? (3) what is the 
correlation between cluster size and thermodynamic stability and is 
there any magic-sized B N cluster with particularly high stability? To 
answer these questions, the low-energy configurations and stability 
of small B N clusters up to N = 20 on Cu(lll) surface were 
systemically investigated by DFT calculations. The ground-state 
structures of Cu-supported B N clusters are shown in Figure 2, and 
some selected isomers are presented in Figure S6 of Supporting 
Information. 

To search the most stable configurations, we have constructed 
numerous structural isomers for each size of B N clusters on Cu sub- 
strate. At least five isomers were considered for N S 9, whereas the 
number of isomers increases with the number of boron atoms, e.g., 
sixteen isomers for B 19 . According to our test calculations, the tri- 
angular networks embedded with quadrilaterals or pentagons (e.g., 
9-b, 11-b, 13-d, 17-e and 19-d in Figure S6 of Supporting Infor- 
mation) have relatively higher formation energies and thus are not 
further considered. Due to the large number of possible conforma- 
tions, we cannot exhaust all isomers for N > 10. Only planar struc- 
tures composed of triangular units with or without hexagonal hole 
were considered, which can be viewed as the precursors of ht-BS and 
f-BS, respectively. In addition, some selected 3D clusters (e.g., ico- 
sahedron for B 12 , double ring for B 20 ) were also placed on Cu( 1 1 1 ) to 
examine the competition between 2D and 3D structures. 

For a single boron atom adsorbed on Cu(lll) surface, the most 
energetically favored position is the bridge site (see Figure 2), whose 
formation energy is 0.09 eV/atom lower than that of the hollow site. 
A boron adatom is unstable on the top site and would move to the 
bridge site upon relaxation. Hence, we tend to avoid the top sites and 
maximize the bridge sites when we constructed the initial geometries 
for Cu(l 1 l)-supported B N clusters. In the case of B 2 dimer, however, 
two boron atoms sit on the hollow sites of Cu(l 1 1) surface with B-B 
distance of 1.60 A since the nearest neighboring bridge sites are too 
close (1.28 A) and the second nearest ones are too far (2.22 A). 

On Cu(lll) surface, B 3 is a regular triangle and B 4 is a rhombus 
composed of two triangles. Chain configurations (3-b and 4-b in 
Figure S6 of Supporting Information) were also considered for B 3 
and B 4 , but are higher in energy by 0.049 eV/atom and 0.095 eV/ 
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Figure 2 | Most stable configurations and formation energies (eV per boron atom, given in the parenthesis) of 2D Bjv clusters (JV = 1 -20) on Cu( 1 1 1 ) 
surface. 



atom, respectively. Hence, chain isomers were not considered for 
larger B N clusters with N > 4. The strong preference of 2D compact 
structures based on triangular units in B N clusters is distinctly dif- 
ferent from Cjv clusters on Ni( 1 1 1 ) surface, which prefer chain con- 
figurations up to N = 1 1 47,48 . This is simply because that boron does 
not adopt sp hybridization like carbon. 

Following the triangular structural motif, the lowest-energy con- 
figurations of Cu(lll)-supported B 5 and B 6 are trapezium and 
parallelogram composed of three and four triangular units, respect- 
ively. The hexagonal ring of B 6 (6-b in Figure S6 of Supporting 
Information) is unstable with higher formation energy by 0.061 
eV/atom, implying that an individual hexagonal boron ring does 
not exist on Cu(lll) surface. An atom-centered hexagon emerges 
as the ground state of B 7 on Cu(lll). Interestingly, the most favored 
configurations for small B N clusters (N £ 7) on Cu substrate are 
all identical to those in vacuum 34 (see Figure S5 of Supporting 
Information for the ground state structures of gas-phase boron 
clusters). 

The ground state structures of gas-phase B 8 and B 9 clusters are 
atom-centered heptagon and octagon (Figure S5), respectively 34 . On 
Cu(lll) surface, however, the wheel structures for B 8 and B 9 are 
unstable and would transform to B 7 -based triangular networks, i.e., 
atom-centered hexagons with one or two atoms on the edge. This 
suggests that the maximum coordination number of boron atom in 
the metal-supported clusters or monolayer sheet cannot exceed six. 
As for B 10 , the lowest-energy geometries in vacuum and on Cu sur- 
face are identical, i.e., a belt-like configuration with regular triangular 
network. Clearly, the most stable structures of small B N clusters (3 £ 
N £ 1 0) on Cu( 111) surface are all made up of triangular units, which 
can be considered as embryo of t-BS. 

On Cu(lll) substrate, B n is the first cluster with triangular net- 
work consisting of a hexagonal hole, which can be regarded as a 
precursor of ht-BS. By contrast, the gas-phase B n configuration is 
constituted by nine triangle units and one quadrilateral 35 . A 2D 
closed-shell structure with C 3V symmetry, consisting of thirteen 



triangles, is the ground state for both gas-phase 34 and Cu-supported 
B 12 . To examine possible existence of 3D structures, B 12 icosahedron 
was placed on Cu(lll) surface (12-c in Figure S6 of Supporting 
Information); but its formation energy is much higher than the 2D 
C 3v configuration by 0.444 eV/atom. Thus, an individual icosahed- 
ron (which is the basic structural unit of boron solids) cannot occur 
on Cu(lll) surface. In previous calculations by Yakobson's group 40 , 
B 12 icosahedron on Ag(lll) surface was also energetically unfavor- 
able with regard to the planar structure by 0.225 eV/atom. 

Similar to B 10 , B 13 on Cu(lll) surface also adopts belt-like con- 
figuration with triangular network. In vacuum, the same structure 
was found for anionic B 13 ~ cluster 34 , but not for neutral B 13 . This can 
be explained by the saturation effect of metal surface since that B 13 
cluster gains about 1.2 electrons from Cu substrate according to 
Bader analysis. The lowest-energy geometry of B 14 is also belt-like 
simply by adding one atom to the B 13 belt (13-a in Figure 2). The 
isomers as pieces of ht-BS were also considered for B 13 and B 14 (13-e 
and 14-c in Figure S6 of Supporting Information), but their energies 
are higher than the belt structures by 0.058 eV/atom and 0.030 eV/ 
atom, respectively. 

In the case of B 15 , ht-BS structural pattern prevails again and its 
lowest-energy geometry is based on 14-c isomer of B 14 with one 
additional atom on the rim, forming a C 2v structure with a hexagonal 
hole symmetrically surrounded by thirteen triangles. Following B 10 
and B 13 , B 16 adopt belt-like triangular network (16-a in Figure 2) as 
its ground state both on Cu surface and in vacuum. Meanwhile, 16-b 
isomer (Figure S6 of Supporting Information) with one hexagonal 
hole is less stable by 0.036 eV/atom probably due to the two-coordi- 
nated boron atom on top of 15-a configuration. 

As a continuation of ht-BS motif, the lowest-energy structure of 
B 17 on Cu(lll) substrate can be obtained by adding a boron dimer 
on top of 15-a configuration, preserving C 2v symmetry. Interestingly, 
it can be regarded as a fragment of every ht-BS considered here 
except the snub one. Similarly, B 18 prefers a triangular network with 
a hexagonal hole in the center and possesses a high symmetry of C 6V . 
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This planar boron patch also exists in most ftf-BSs except g 2/15 -sheet 
and pVsheet. The 18-c isomer as a full triangular network with D 3h 
symmetry, which is the ground state in vacuum, is higher in energy 
by 0.027 eV/atom on Cu surface. At N = 19, the triangle-based belt 
configuration prevails again and its energy is lower than the 19-b 
isomer of ht-BS motif (Figure S6 of Supporting Information) by 
0.038 eV. 

The most stable structure of B 20 on Cu surface can be obtained by 
adding two boron atoms to the 18-a configuration, resulting in C 2y 
symmetry. Again, it can be viewed as a fragment of g 1/8 -sheet, a- 
sheet, al-sheet, and struc-l/8-sheet. For comparison, the 3D double- 
ring tubular structure (20-e in Figure S6 in Supporting Information), 
which is the ground state of B 2 o in vacuum 33 , is energetically unfa- 
vorable on Cu surface with AE = 0.073 eV/atom, simply because 
only half of boron atoms interact with Cu atoms. The high formation 
energies of B 12 icosahedron and B 20 double-ring clearly indicate that 
3D structures are not favored on Cu(lll) surface. 

From the above results, we can derive some general trends about 
the structures and growth patterns for small boron clusters on 
Cu(lll) surface. First, all Cu-supported B N clusters prefer 2D con- 
figurations rather than 3D ones. Such 2D preference is certainly a 
prerequisite for fabricating boron monolayer films on metal sub- 
strates. At cluster sizes of N = 3, 4, 5, 6, 7, 10, 12, 16, the lowest- 
energy geometries of B N in vacuum (highlighted in Figure S5 of 
Supporting Information) and on Cu surface are identical and they 
are all symmetric triangular networks. However, this kind of geo- 
metries cannot form in the B u , B 13 , B 15 , B 17 and B 19 clusters; a 
quadrilateral has to be included in the gas-phase configuration to 
reduce the number of low-coordinated boron atoms. When sup- 
ported on Cu( 111) substrate, a hexagonal hole instead of a quadrilat- 
eral is formed inside the triangular networks of B N clusters. These 
structures can be considered as the early stage of ht-BS. Therefore, we 
infer that metal surface plays a significant role in occurrence of the 
boron sheets composed of triangular units and hexagonal holes (i.e., 
ht-BS). 

For B N clusters with N s 11, we observed strong competition 
between f-BS and ht-BS motifs. For those full triangular islands as 
precursors of f-BS, belt-like structures are more favorable and 
emerge atN= 10, 13, 14, 16, 19, whereas only B 12 prefers an enclosed 
triangular network. On the other hand, in those lowest-energy tri- 
angular networks with hexagon hole (11-a, 15-a, 17-a, 18-a and 20-a 
in Figure 2), the hexagons are always well embraced by triangular 
units. We further constructed several belt-like configurations incor- 
porated with one hexagonal hole (12-b, 15-b and 18-b in Figure S6 of 
Supporting Information) and 2D triangular networks with two hexa- 
gonal holes (14-d and 16-e in Figure S6 of Supporting Information). 
All of them were found to be energetically unfavorable. Therefore, we 
can conclude that a hexagonal hole has to be enclosed by sufficient 
amount of triangles to achieve stable 2D boron cluster within ht-BS 
motif. 

Interaction between boron clusters and Cu(lll) surface. We first 
discuss interaction between 2D boron clusters and Cu(lll) surface 
in terms of buckled heights (Figure S7 of Supporting Information). 
For three typical enclosed triangular boron clusters (12-a, 16-c, 19-e) 
as precursors of f-BS, the buckled height increases with increasing 
cluster size and reaches up to 1.32 A at N = 19, compared to h = 
1.21 A for the infinite triangular sheet. The buckled heights for two 
representative hole-doped triangular boron clusters (17-a and 18-a) 
oi ht-BS type are only 0.43 A and 0.28 A, respectively, much less than 
that of the infinite ht-BS (about 1 A). The larger buckling in Cu- 
supported ?zf-BS is mainly caused by the mismatch of 2D boron 
lattice and Cu(lll) surface. In contrast, small boron clusters with- 
out constraint of periodic boundary condition can accommodate 
Cu(lll) surface better since the edge boron atoms have more 
degrees of freedom for adjusting their adsorption positions. 



As illustrated by charge density difference (Figure 3a) and Bader 
charge analysis (Figure 3b), interactions between boron clusters and 
Cu( 111) surface mainly occur on the periphery of boron clusters. Cu 
substrate causes depletion of in-plane ct states and accumulation of 
out-of-plane n electrons on the borders of boron clusters. The 
amount of electron transfer from Cu substrate to boron clusters is 
usually larger than that for boron monolayers (about — 0.04 electrons 
per B atom). The larger electron transfer in Cu-supported Bjvdusters 
remedies the electron deficiency of boron and thus results in different 
equilibrium geometries from the gas-phase clusters. 

Previously, the extraordinary stability of free-standing a-sheet of 
boron were explained by the 3c-2e a-bonds, 4c-2e a-bonds, and 
delocalized n -bonds (6c-2e) 20 . Due to lack of periodic environment 
and the consequently unsaturated periphery atoms, B N clusters as 
small fragments of ht-BS, e.g., 17-a and 18-a, are not energetically 
favorable in vacuum but can be stabilized on Cu(lll) surface. This 
can be understood by the following bonding picture. There are one 
3c-2e CT-bond in each triangle and one delocalized 7t-bond in each 
hexagonal hole, respectively. On the other rim of supported boron 
clusters, the corner boron atoms form two bonds with metal sub- 
strate, while the other edge atoms form one bond. For example, B 18 of 
18-a configuration has eighteen triangles and one hexagon, which 
require 18 X 2 = 36 <r electrons and 2 n electrons. On the periphery 
of B 18 , each of the six corner B atoms forms two bonds with Cu 
atoms, and each of the six edge B atoms forms one bond, respectively. 
Hence, totally 56 electrons are needed to fulfill all these bonds. 
Indeed, eighteen boron atoms can offer 54 electrons and the entire 
B ls cluster gains about 1.8 electrons from Cu surface according to 
Bader analysis, making the total number of electrons near 56. Clearly, 
passivation effect by metal surface and charge transfer from metal to 
boron account for stabilization of the ht-BS patches. 

Implication on the growth behavior of boron monolayer. Figure 4 
plots formation energies of free-standing and Cu-supported B N 
clusters as function of size N. Due to metal passivation, the 
formation energies of B N clusters are significantly reduced from 
vacuum to Cu(lll) surface, especially for those small clusters. For 
example, the difference of formation energies for systems before and 
after adsorbed on Cu(lll) surfaces is 1.359 eV/atom for B 6 and 0.5 
~ 0.7 eV/atom for N = 12 ~ 20, respectively. As clusters further 
increases, it eventually reaches about 0.15 ~ 0.19 eV/atom for 
infinite /if-BSs (see Table 1 and Figure 4), e.g., 0.188 eV/atom for 
g 2/15 -sheet. Again, this indicates that large boron sheets interact with 
Cu substrate more weakly than small B N clusters, in line with above 
charge analysis. It is noteworthy that the cohesion strength between 
various boron sheets and Cu(lll) surface is the same magnitude of 
graphene on Cu surfaces (about 0.022 eV/A 2 ) 50 , implying that the as- 
fabricated ht-BSs may be peeled off from metal substrates like 
graphene. 

Within the explored size range, the formation energy of B N clus- 
ters decreases as N increases both in vacuum and on Cu(l 1 1) surface, 
excepted for some highly stable sizes. Supposing boron atoms are 
directly deposited on Cu( 111) surface, the growth process of a mono- 
layer boron sheet can thus be viewed as a phase transition from 
dispersed boron monomers to 2D sheet with patterned structures. 
Therefore, the chemical potential of various BSs on Cu(lll) at 0 K 
can be computed by: 

A/i = Ef orm (B 1 ) —Ef orm (BS), (1) 

where Ef orm (B!) = 1.993 eV/atom is the formation energy of boron 
monomer on Cu(lll) surface, Ef orm (BS) is the formation energy of 
boron sheet on Cu(lll) surface, which lies between 0.207 eV/atom 
and 0.338 eV/atom. Thus the value of Au changes from 1.786 eV/ 
atom to 1.655 eV/atom for various patterned structures of BSs on 
Cu(lll) surface. 
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Figure 3 | (a) Charge density difference and (b) on-site charge transfer fromBader analysis forfour typical 2DB N clusters (N= 12, 16, 17, 18) onCu(lll) 
surface. In (a), electron accumulation is represented by green and electron depletion is represented by blue. In (b), the unit of charge transfer is lei. 



Using chemical potential defined in Eq.(l), the Gibbs free energy 
of boron clusters on Cu(lll) surface can be calculated 51 by 



AG = AE B -NxAfi, 



JE B =Nx E form (B N )-Nx E form (BS) 



(2) 



(3) 



where N is the number of boron atoms in cluster, E form (B N ) is the 
formation energy of boron cluster, Ef orm (BS) is the formation energy 
of infinite boron sheet. As displayed in Figure 4b, the Gibbs free 
energy B N cluster is always negative and decreases monotonously 
with increasing cluster size N, which implies spontaneous aggrega- 
tion of boron atoms and incessant growth of boron cluster on 
Cu(lll) surface at normal conditions. In comparison, on Ag(lll) 
surface, AG slowly decreases with increasing cluster size even for 
small A(i = 0.2 eV/atom and the nucleation barrier is about 3 eV 40 . 

Besides, some B N clusters (e.g., B 7 , B 10 , B 12 and B 16 ) exhibit high 
stability on Cu substrate from the second derivatives of the formation 
energies (see Figure S8 of Supporting Information), and they possess 
the same planar structures as in vacuum. Therefore, these boron 2D 
clusters can be directly soft-landed on Cu(l 1 1) surface using cluster 
beam apparatus 38 as nucleation seeds for further growth of various 
boron sheets. 

In principle, boron sheets on Cu(lll) surfaces can be obtained 
via continuous growth of small clusters or coalescence of them. 
According to our calculations (Figure S9 of Supporting Infor- 
mation) , the diffusion barrier of isolated B atoms on Cu( 1 1 1 ) is only 
0.141 eV. Therefore, the continuous growth of 2D boron clusters is 
guaranteed by the negative AG and low diffusion barrier. To schem- 
atically show how a ht-BS forms on metal surface, in Figure 5 we 
propose two possible growth pathways based on lowest-energy or 
metastable configurations of B N clusters. In the path (a), starting 
from N = 11, small B N clusters with hexagonal-hole-doped triangu- 
lar lattice occur on Cu surface as nucleation centers of ht-BS. As 
boron atoms aggregate around the B N cluster, the hexagonal hole 
is completely surrounded by triangular units at B 18 . Afterwards, 
more boron atoms will accumulate around B 18 , e.g., a B 2 dimer on 
top of B 18 leading to B 20 - Although larger B N clusters (N > 20) are not 
explored due to restriction of computing resource, one can naturally 
anticipate formation of more and more hexagons inside the triangu- 
lar network of a 2D B N cluster as it grow bigger, finally coming into 
being macroscopic monolayer of ht-BS. 



As shown in the path (b) of Figure 5, boron nanoribbons may also 
serve as building blocks of ht-BS. Since the belt-like B N clusters (i.e., 
Bio. B 13 , B 16 , B 19 ) with regular triangular lattices are very stable, 
many nanoribbons of finite length may exist during the initial stage 
of growth. When several boron ribbons aligning in different orienta- 
tions merge with each other, a hexagonal hole may occur via either 
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Figure 4 | (a) Formation energies of the most stable B N clusters and boron 
sheet in vacuum and on Cu( 1 1 1 ) surface, (b) The Gibbs free energy of Cu 
supported boron clusters as a function of cluster size. 
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Figure 5 | Schematic plots for two possible growth paths of hf-BS (g 2 /i5-sheet as representative) from small B^clusters. Formation energies in unit of eV 
per boron atom are given below the name of boron clusters and sheet. 



directly enclosing of three nanoribbons with well-defined intersec- 
tion angle of 120° or diffusion of few boron atoms (see discussion 
below for the formation mechanism of hexagonal hole). 

Is it possible to grow a boron monolayer fully with triangular 
lattice on Cu substrate? B 12 as a possible precursor of t-BS (12-a in 
Figure 2) is rather stable on Cu(lll) surface. However, as the cluster 
becomes bigger, these enclosed triangular structures resembling B 12 
(14-b, 16-cand 19-e in Figure S6 of Supporting Information) become 
less favorable. The 19-e isomer with a triangular network of C 6V 
symmetry forms a dome-like shape on Cu surface similar to gra- 
phene patch 49,52 , and its formation energy is 0.121 eV/atom higher 
than the 19-a ground state. The higher formation energies of 
enclosed triangular boron clusters originate from their larger buckled 
height h, which are shown in Figure S7 of Supporting Information 
(B 12 , B 16 and B 19 ), and prevent formation of triangular boron sheet. 

To achieve the patterned fof-BSs (in Figure 1 and Figure SI of 
Supporting Information) from continuous growth and coalescence 



of small boron clusters, creation and redistribution of hexagonal 
holes are the crucial steps, which can be experimentally realized via 
proper annealing. To explicitly see how a hexagonal hole arise within 
a finite triangular network and further diffuse inside a B N cluster, we 
performed NEB search for the structural transformation from tri- 
angular network to hexagonal hole-doped triangular network using 
B 15 as an example. As shown in Figure 6, the initial configuration of 
B 15 is obtained by adding one additional atom to the metastable 
structure of B 14 (14-b) on Cu surface; but this fully triangular net- 
work is not favorable for B 15 . After one atom (highlighted in Figure 6) 
jumps to the adjacent triangular site, a hexagonal hole occurs as an 
intermediate (IM) state. The activation barrier for creating of a hexa- 
gonal hole as the transition state one (TS1) is 0.531 eV. When a 
hexagonal hole arises at the edge of a triangular network, it can 
further diffuse into the cluster interior to be enclosed by as more 
triangles as possible. As depicted in Figure 6, after the highlighted 
boron atom migrates from the edge of hexagonal hole to its center 




Ini: 0.000 eV 




Fin: -0.832 eV 



Figure 6 | Schematic plots for formation and diffusion of a hexagonal hole inside B 15 cluster. The energies for the initial (Ini), transition (TS1, TS2), 
intermediate (IM), and final (Fin) states are given. The migrated atom and its diffusion direction are highlighted. 
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with a small activation barrier of 0.49 1 eV passing through the trans- 
ition state two (TS2), the energy of Cu-supported B 15 cluster is 
reduced by 0.832 eV. Such small energy barriers for creating and 
diffusion of a hexagonal hole ensure that a hexagonal hole can arise 
simultaneously once the triangular network is energetically unfavor- 
able. In principle, the driving force for formation/diffusion of hexa- 
gonal holes inside a triangular patch during growth process 
originates from the formation energy differences AEf orm of various 
boron sheets, i.e., AEf orm 2 0.1 eV/atom between /if-BSs and (-BS, 
AEform = 0.039 eV/atom between g 2 /i 5 -sheet and struc- 1/8-sheet 
(see Table 1). 

Discussion 

Various types of boron monolayers on Cu(lll) surface were inves- 
tigated using comprehensive ab initio calculations. Most boron 
monolayers are very stable on Cu(lll) surface and retain the same 
geometries as in vacuum. The formation energies for these boron 
monolayers are reduced on Cu surface due to metal passivation 
effect. The final boron sheet synthesised on Cu surface would be 
the mixture of these /jf-BSs due to the small energy differences 
between these 2D isomers. AIMD simulation at 500 K and 800 K 
comfirms high thermal stability of the boron sheets on Cu(lll) 
surface. 

The lowest-energy and metastable configurations of boron clus- 
ters B N (N £ 20) on Cu( 111) surface were systematically investigated 
and compared with those in vacuum. On Cu(lll) surface, boron 
atoms tend to form 2D compact clusters on the basis of triangular 
units and their formation energies continuously reduce with cluster 
size, which imply the possible incessant growth of 2D boron cluster 
on Cu substrate. There is a competition between fully triangular 
structure and mixed triangular and hexagonal structure. Some frag- 
ments of ht-BS can be very stable due to metal passivation effect and 
charge transfer from metal surface to boron clusters, which offer the 
prerequisite for synthesis of ht-BS. The small diffusion barrier 
(0.141 eV) also suggests fast migration of boron atoms on Cu surface, 
which is beneficial for growth of 2D boron patches. NEB search 
shows that hexagonal holes can easily arise and diffuse in boron 
cluster during the growth on Cu(lll) surface with relatively low 
energy barriers (about 0.5 eV). Therefore, we proposed two possible 
ways of boron monolayer growth, i.e., either depositing boron atoms 
directly on Cu( 111) surface or soft landing of small planar B N clusters 
generated by cluster beam apparatus. Both approaches result in 
mixed hexagonal-triangular boron sheet on Cu(lll) surface after 
nucleation and thermal annealing. In order to obtain boron mono- 
layer, the evaporation rate of the boron source and the total amount 
of boron atoms deposited on metal substrate must be strictly con- 
trolled. The present theoretical results are undoubtedly very helpful 
for understanding the atomic structures and growth mechanism of 
boron monolayer on Cu(lll) surface and provide a useful guideline 
for experiment synthesis of monolayer boron films. 

Methods 

All calculations were carried out using the Vienna Ab initio Simulation Package 
(VASP) based on density functional theory (DFT) and planewave basis 53 . The elec- 
tron-ion interactions were described by the projector augmented wave (PAW) 
potentials 54 . To treat the exchange- correlation interaction of electrons, we chose the 
Perdew- Burke- Ernzerhof (PBE) functional within the generalized-gradient 
approximation (GGA) 55 . A kinetic energy cutoff of 400 eV for the planewave basis 
and a convergence criterion of 10" 4 eV for the total energies were carefully tested and 
adopted for all DFT calculations. 

Here we chose Cu(l 1 1) surface as substrate since that it is easy to build co-periodic 
lattice of boron monolayer and Cu(lll) surface and that the boron atom has mod- 
erate adhesion energy on Cu surface (about 0. 1 eV/A 2 ) 40 . Meanwhile, Cu( 111) surface 
is a widely used substrate for epitaxial growth of graphene 56,57 . The Cu(lll) surface 
was modeled by a three-layer slab model within 2D periodic boundary condition, 
which was cleaved from fee Cu solid with an equilibrium lattice constant of 3.63 A 
from our first-principles calculations. To avoid interactions between adjacent peri- 
odic images, a vacuum space of more than 12 A was included in the slab model. With 
fixed supercell parameters, the three-layer slab model was further relaxed with the 



bottom layer fixed to mimic a semi-infinite solid. The validity of this three-layer slab 
model was assessed by our test calculations with a five-layer slab model. 

Following recent progresses on the 2D boron sheets 1316 ' 17 ' 21-24 , we considered nine 
possible configurations of boron monolayer on Cu{lll) surface. The co-periodic 
lattices of these boron monolayers and Cu(l 1 1) surface were built by compressing or 
stretching the Cu lattice slightly to fit the 2D unit cell of boron sheet. The small lattice 
mismatches (A a and Ab along two crystallographic orientations) between boron 
sheets and substrates are denned as: 

d a — \dsheet ~ ^sub \ , (4) 

Gsub 



db = I \bskeet — b su y\, (5) 

Osub 

where a S hect> bsheet* a sub ana< hub are tne supercell parameters of the boron sheet and 
the Cu substrate, respectively. The detailed structural information for these co-peri- 
odic lattices of boron sheets and Cu substrate are summarized in Table Si of the 
supporting information (Supporting Information). Note that current supercell 
models of co-periodic lattices are a compromise of computational cost and lattice 
mismatch. In reality, there should be almost no mismatch if the co-periodic supercell 
is big enough. 

In addition to the infinite boron sheets, B N clusters of different sizes (up to N = 
20) were placed on the slab models of Cu(lll) surface and the entire cluster- 
substrate systems were fully relaxed (also with the bottom layer of Cu atoms 
fixed). The lateral dimensions of the supercell for these cluster- substrate systems 
were chosen to be sufficiently big to ensure the distance between adjacent clusters 
no less than 10 A. 

During geometry optimization, the spacing of k point grids was chosen to be 

0. 03 A" 1 and a convergence criterion of 0.02 eV/A for force was adopted. In order to 
obtain more accurate energies and electronic properties, denser k-point meshes with 
uniform spacing of 0.02 A" 1 were used for calculations of total energy, charge density 
difference and charge transfer. Bader analysis 58-60 was preformed to evaluate the 
charge transfer between Cu(l 1 1) surface and boron sheets or clusters. 

The strength of interaction between boron monolayer and Cu( 111) surface can be 
characterized by the cohesive energy E c defined as: 

E c = — (E su b + E sheet —Ef), (6) 

where E su y is the energy of substrate, E s h ect is the energy of boron monolayer, E t is the 
total energy of boron monolayer on Cu(lll) system, N is the number of boron atoms. 

To describe the stability of a boron monolayer or cluster on Cu( 111) surface, we 
define its formation energy as: 

E form = - (Et -E sub -N x £ B ), (7) 

where E t and E su (, are the same as the definitions inEq. (6),E B is the energy per atom in 
the boron solid of a phase 44 . NTs the number of boron atoms in boron monolayer or 
B w cluster. In a similar manner, we can define the formation energy of boron 
monolayer or cluster in vacuum {Ef orm vac ) using Eq. (7), where E t becomes the total 
energy for boron monolayer or cluster only, and the E su b term vanishes. 
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